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Abstract. We discuss atomic lattice excitons (ALEs), bound particle-hole pairs 
formed by fermionic atoms in two bands of an optical lattice. Such a system 
provides a clean setup to study fundamental properties of excitons, ranging from 
condensation to exciton crystals (which appear for a large effective mass ratio 
between particles and holes). Using both mean- field treatments and ID numerical 
computation, we discuss the properities of ALEs under varying conditions, and 
discuss in particular their preparation and measurement. 



PACS numbers: 03.75.Ss, 71.35.-y, 42.50.-p 
1. Introduction 

Excitons, the bound state of a fermionie particle and a fermionic hole, belong to 
the most basic excitations in any system of fermionic particles confined in the band 
structure of a periodic potential [HE]- In a semiconductor, metastable excitons are 
created by optically pumping electrons from the valence to the conduction band, and 
Coulomb repulsion between the electrons then leads to the formation of a bound 
state of the electron with the hole it left behind. At sufficiently low temperatures, 
the approximately bosonic statistics of the composite objects give rise to a superfluid 
state corresponding to a condensate of excitons [TJ[51[31[31[5J[ni[71[Hl[ni HH] ■ 

In the present work we investigate the properties of bound particle-hole pairs 
of fermionic atoms confined in an optical lattice, where two-particle interactions are 
provided by collisional processes between atoms [TTl [12l [131 Hi] > or Atomic Lattice 
Excitons (ALEs) [H]. Optical Lattice systems [HIISlinillHlIISlEniEIlEaES], 
have emerged as a very attractive tool for the study of many-body lattice models 
(see, e.g., O [25l [Ml Ell EH EH [30] particularly because of the high degree 
of available control over system parameters (e.g., control over interactions using 
Feshbach resonances [31] |32l [33l [34]. These systems also allow for a wide 
range of measurement techniques, giving access to quantities including momentum 
and quasimomentum distribution, and higher-order correlation functions via noise 
correlation measurements [35^ |36l [37] . 
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Figure 1. Atomic Lattice Excitons (ALEs): (a) ALEs can be prepared from a 
fermionic band insulator, with one atom per site in the valence band (denoted 
V). (b) Pairs of conduction band atoms and valence band holes are created by 
transferring valence band atoms to the conduction band (denoted C), which is 
separated from the valence band by the band gap energy Eg. c) Particle- hole 
pairs move by tunnelling of the conduction-band atom with rate Jc and back 
tunnelling of a valence-band atom with rate J„ , which corresponds to co-tunnelling 
of the hole. Co-tunnelling occurs because of the interband atom-atom repulsion 
U. Thus, atom-hole pairs move together, forming the exciton. 



ALEs constitute a novel realisation of excitons, allowing fundamental aspects 
of excitons as composite objects and as an interacting many body system to be 
investigated in a particularly clean and controlled setting. In fact, the many body 
physics represented by ALEs corresponds closely to the basic theoretical models, 
which have been developed in seminal work in the context of exciton condensates 
in semiconductors physics, focusing mainly on many body physics of excitons 
as interacting composite objects of two fermionic electron-hole constituents. A 
semiconductor environment constitutes, of course, a much more complex system with 
coupling to other degrees of freedom, e.g., lattice phonons and vacuum modes of 
the radiation field, which are essentially absent in atomic lattices. Furthermore, 
additional parameter regimes can be accessed for ALEs, producing quantum phases 
that are not normally observed in semiconductor systems. This includes, for example, 
the formation of exciton crystalline structures in regimes where the effective mass of 
particles and holes on the lattice are substantially different. 

In the simplest case, ALEs can be formed in a system of spin-polarised fermions 
in a optical lattice, where the lowest Bloch band (called the valence band) is initially 
fully occupied (figured^). This can be set up in ID, 2D or 3D, although in this article 
we will primarily concentrate on the ID case. This is both the conceptually simplest 
case as only a single excited Bloch band needs to be considered, and also provides 
interesting analogies to the possible formation of excitons in quantum wires. The 
lattice is assumed to be sufficiently deep for a tight-binding approximation I16j to be 
valid for the first two Bloch bands, which are separated by the band gap Eg. Atomic 
lattice excitons (ALEs) are created from the ground state by exciting atoms with a 
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laser driven Raman process into the next highest Bloch band (called the conduction 
band), where they have a tunnelling rate Jc- The excitation leaves a fermionic hole 
in the valence band, (figure [T|d), moving with hopping rate J/, = — J^, where is 
the valence band hopping rate of the atoms. While spin-polarised fermionic atoms 
have, by symmetry, no s-wave scattering, a repulsive coUisional interaction U can be 
generated between conduction and valence-band atoms in the same lattice site (e.g., 
by using an off-resonant Raman transition to mix in some component of a different 
spin state for particles in one of the two bands.). This on-site repulsion, in turn, gives 
rise to an effective attraction between conduction-band atoms and holes. 

In this work we investigate ALEs, applying both mean- field theory and exact 
numerical calculations in ID. We derive the Fermi-Hubbard Model description of 
ALEs, and then discuss both single ALEs (Sec. Wh) and interaction between two 
ALEs (Sec. [2j3). In particular, interactions between long-lived ALEs, which due to 
the large effective mass ratio between particles and holes can be analysed in a Born- 
Oppenheimer approximation, are characterised by a finite range repulsion. In Sec. [3] 
we then discuss the resulting many-body phases, including ALE condensates, and a 
crystal phase that arises from the effective long-range interactions. We discuss how 
the properties of these phases can be measured using RF Spectroscopy and noise 
correlation techniques (Sec. H]), and discuss schemes to prepare well-defined filling 
factors of excitons in a low energy many-body state in an optical lattice (Sec. [5]). We 
give a summary in Sec. [B] In [Appendix A we provide extra details of our solution for a 
single exciton on the lattice, while in Appendix B| we detail the solution to the problem 
of A'^ static lattice impurities that underlie our Born-Oppenheimer approximation. 
Finally, in [Appendix C| we provide additional information about the general treatment 
of excition condensation. 



2. Atomic lattice excitons (ALEs) 

We consider a gas of ALEs, formed in an optical lattice that is sufRciently deep for a 
tight-binding approximation [16] to be valid for the first two Bloch bands, which are 
separated by the band gap Eg. If the lattice is isotropic or near- isotropic, there will 
be degenerate lowest excited p- bands, as discussed in [27]. In this paper we focus on 
excitons that are tightly confined in two dimensions so that we consider excitons in ID, 
and will also give mean-field results that are applicable to excitons in 2D or 3D with 
strong anisotropics so that the p-band degeneracy is lifted, and one of the bands can be 
chosen as the conduction band. In the conduction band, atoms have a tunnelling rate 
between neighbouring sites with lattice vectors x and x' (these tunnelling rates 
can, in general, be anisotropic). The excitation then leaves a fermionic hole in the 
valence band, (figure [Dd), moving with hopping rate J]^^ = —J^^ , where is the 
valence band hopping rate of the atoms between neighbouring sites. It is convenient 
to expand the field operators of the lattice-confined fermions in the localised Wannier 
modes w'f'^'''^(r — x) and the corresponding annihilation operators Cx'^'"^ of conduction 
and valence band respectively [16] 

^{c..}(r)^^4c,.}^{c,.}(r_x), (1) 

X 

where x = a(xi, . . . tXd), (a;i, . . . tXd) G 2^ and a is the lattice spacing. 

A repulsive interaction U between conduction and valence-band atoms in the 
same lattice site can be generated as discussed in sec. [TJ which then gives rise to 
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an effective attraction between conduction-band atoms and holes. A Hamiltonian 
describing the interacting populations of atoms in two bands can be obtained by 
inserting the decomposition[l]into the two-species Hamiltonian in second quantization. 
We obtain a repulsive two-species Hubbard-Hamiltonian, 

^ = - E •^c'^'^'^cx' - E J^^'^W + uY.^ibib^^^^ (2) 

(XX') (yy) X 

where cl^ and are atom creation operators on lattice site x in conduction and 
valence-band respectively, obetying the usual fermionic commutation relations. This 
Hamiltonian is valid in the limit where U, Jc, Jy Eg, where Eg is the band gap. This 
is well fulfilled by typical experimental parameters. Note that terms involving transfer 
of particles from one band to another do not appear because of this condition, and 
therefore the total number of valence band atoms and the total number of conduction 
band atoms are conserved separately. 

Holes appear through the introduction of fermionic hole creation operators d^. 
in the lower band by the definition d^. = b^. This trivial substitution yields the 
particle-hole Hamiltonian 

^ = - E Jp'^'^l^^' - E Jh'^'dldy -UY, 44d.Cx (3) 
{xx') (yy) X 

where — , — — J^""" are now the hopping rates of conduction band 

particles and holes respectively, and U is the strength of the particle-hole attraction, 
for conduction band atoms and valence band holes, respectively The attractive 
interaction may give rise to bound states of conduction band atoms and valence band 
holes, i.e., ALEs. Although we do not consider this case here, note that U < would 
give rise to repulsively bound excitons, in the same sense as repusively bound atom 
pairs, which were recently observed experimentally |17| (These would consist of bound 
pairs where the bound state is higher in energy than the continuum of scattering states, 
instead of the stable bound states where the bound state has a lower energy than the 
scattering continuum). In the case of excitons confined to move in ID, we use the 
notation Jp — J^-^ and Jh — J^'^ ■ In the setup we consider, we will have Jh < 
and Jp < 0. 

2.1. The single ALE 

In order to obtain the wavefunction for a single ALE, we make an ansatz for the 
exciton creation operator: 

itb)-E<^-y44l«)' (4) 
xy 

where \v) is the particlc-holc vacuum, and 4>x.y is determined from solution of the 
eigenvalue equation HA^\v) — EA^v). We will now investigate the form of cjix.y 
for the case of excitons confined in ID, a result which can also be straightforwardly 
generalised to 2D or 3D [T71[3H]- 

We can obtain a simple analytical solution that has the same form either for the 
case Jp = Jh, or the case of strongly imbalanced hopping rates, |Jp| 3> |J;i|. In the 
first case, introducing Center-of-mass (COM) and relative coordinates R = {x + y)/2, 
r = X — y gives rise to the product ansatz 

(l)x,y = exp [iK{x + y)/2] p^^y, (5) 
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which solves the eigenvalue equation exactly. The COM quasimomentum K e 
[— 7r/a, TT / a] is factored out, and we are left with equations for the relative wavefunction 

Pr- 

~ J{Pr+l + Pr-l) ~ USrOPr = Epr- (6) 

Here J = Jk — "^Jp cos{Ka/2) is the effective hopping rate in relative coordinates and 
6ro the Kronecker-delta. Note that the relative wavefunction pr still depends on the 
centre-of-mass quasimomentum K. 

The parameters encountered in an experiment for the suggested implementations 
of ALEs will typically fulfill |^?!,|/|^p| ^ 1, so that holes will appear much heavier 
than the conduction band atoms. We can make the approximation that the centre of 
mass is approximately located at the position of the hole, and obtain in this sense a 
Born-Oppenheimer wavefunction for the single exciton, 

<Px,y = C{y)p^-y, (7) 
where C{y) is the wavefunction for the hole. This leads to[6l but now with J = Jp. 

Eq. [6] can be solved exactly either by direct solution of the difference equation 
in ID using a standard exponential ansatz, or using Green's function methods 
(see I Appendix A[ ). For the case of Jp = Jh, It yields a single bound state solution for 
each value of the COM quasimomentum and a continuum of solutions describing 
unbound states, illustrated in figure [2^. The continuum solutions correspond to 
scattering states, whilst the bound state solution appears as a Bloch band for the 
composite object, i.e., the ALE. The wavefunction pr for the bound state solution is 
given by 

^^zE^VW+MEEM (8) 

E^-^lPTj2jf (9) 

where C is the normalisation constant and E < the energy of the bound state 
solution. This is depicted in figure Od. As U/\J\ increases, the ALE becomes more 
tightly bound, and p,. decays more rapidly. For J > 0, the phase of the wavefunction is 
constant (given by the phase of C), whereas for J < 0, the sign of the pr will oscillate 
between neighbouring sites. 



2.2. Tightly Bound ALEs 

In the tightly bound limit, \Jp\,\Jh\ ^ U , we can interpret ALEs as effective hard-core 
bosons with creation operators 6j, — cj,dj, [39j . The effective Hamiltonian is given by 

Hhc = -Jcs ^ i>\bj + Foff X! (10) 

{ij) iv) 

where ni — blbj. The effective exciton hopping rate JcS can be calculated in 
degenerate second-order perturbation theory [39j . and the dimension- independent 
result is Jcs — 2JpJh/U. A weak nearest-neighbour repulsion also appears in 
second-order perturbation theory, with Vcs — {Jp -\- JD/U again independent of 
dimensionality [13 SO]- Note that in this limit, where the atoms are treated as hard- 
core bosons, this form for the interaction precludes the formation of bound states (i.e., 
bi-excitons). 

X Note that whilst bound states exist in ID for any positive value of U , if these equations are 
solved in 3D, then a bound state only appears above a critical value for the interaction strength 
{U/\J)crU ~ 3.95. 



Atomic lattice excitons: from condensates to crystals 



6 



CM 



scattering 
xontinuum 



^0 
-U- 



bound particle-hole pair: 
exciton 





Kaj-K 




Figure 2. (a) Energy eigenvalues for a single particle and a single hole on a 
ID lattice plotted as a function of COM quasimomentum, as found from solution 
of l(6ll, with J = (see text). The lower curve represents the Bloch Band 
of a single bound atom-hole pair, the ALE. The upper portion of the spectrum 
is the scattering continuum of the unbound atom and hole, where the shading 
corresponds to the density of states (darker shading for higher density of states), 
(b) Bound-state relative wavefunctions at different binding energies -Efc/J, valid 
both for Jj, = Jp (J = Jjf, see text), and also in the limit <C = J- 
Weakly bound states have an appreciable range over 10 or 20 sites. This falls 
substantially for moderately bound states (dashed line). 



2.3. Interaction of ALE s: The Born- Oppenheimer approximation 



In the typical experimental situation, where Jp ^ Jh, we can treat ALE- ALE 
interactions in a Born- Oppenheimer (BO) approximation, and observe longer range 
interactions mediated by the faster tunnelling conduction band atoms. We find it 
convenient to derive the results in this limit for the general case in D dimensions, 
although we will again primarily apply them in the ID case. 

The basic idea is that we initially assume the holes to be essentially static (i.e. 
Jh ~ 0) at a given set of positions R — {Ri,R2, . . . jRat}. The conduction band 
atoms then move with hopping rate Jp in the static potential given by the holes, and 
have coordinates denoted r = {ri, r2, . . . , r^}- In this sense, we begin by decomposing 
the full Hamiltonian into two parts, one describing slow- moving holes {Hh) and one 
particle motion and particle- hole interaction Hbo- In first quantisation we obtain 



H =Hh 



H 



BO 



N 



n=l 



N 

E 



N 



-JrA 



p^r„ 



r„Rv 



Hbo 

where the operator 

D 

Ax*(x,y) = ^ [*(x+ed,y) + 1'(x-e<j,y) - 2^'(x,y)] 



(11) 
(12) 

(13) 



denotes a discrete lattice Laplacian on a cubic lattice with unit vectors in D 
dimensions. Note that we have shifted the zero of energy by N{2Jp + 2Jh) for 
convenience in writing the discrete Laplacian. 

We write the full time-dependent many-body wavefunction of the ALEs as 



V'(r,R, i) - 51 ^"(^' ^)'^°(£'^)- 



(14) 
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Here, the functions 4>a (z; R) are solutions of the Schrodinger equation for motion of 
the particles in the static potential provided by the holes, 

iJso [E]'/'a(r;R) = ^a;Bo[R]</'a(r; R). (15) 

where R are considered parameters. 

In order to obtain an equation for the functions Ca (R, t) , we apply the full 
Hamiltonian to (HH), multiply with (/)^(r;R) and trace over the rs: 



i9tC„(R,i) 



N 



-Jh ^ ^R„ + -£'q;Bo[R] 



n=l 

N D 



Cq.(R, t) 



-•^'^ E E E E (e. + .t)\Y.^l (i; R) [^P (r; R + - (i; R)](|6) 

p n=lti=le=± [_ r 

where e^;" = {0, . . . , 0, e^;, 0, . . . , 0} represent the D different unit vectors on the 
lattice for the n-th hole coordinate. This describes the dynamics of our system as 
a multi-channel problem, with each a denoting a particular channel, corresponding 
to an eigenstate found from (fT5|l . The first term describes the hole dynamics for a 
single channel, with the holes moving in the effective potential i?Q;_Bo[R]- The second 
term provides coupling between the channels, and has the general form of another 
discrete Laplacian acting on Cp (R, t) , where each hopping process now carries R- and 
a-dependent coefficients. 

If the terms corresponding to coupling between the channels are small, then we 
can reduce the wavefunction expansion, (|14p to a single value of a. This we will do for 
the lowest energy solution to (fT5| . 07(r;R), which corresponds to N bound ALES. 
We can then interpret the corresponding function C'-y (R, t) as the wavefunctions of 
the composite ALEs, moving in a potential obtained from the wavefunctions (^-v(r: R). 
This interpretation amounts to an adiabatic approximation in the ratio of the coupling 
terms in (flBj) with 77^/? and the energy separation |£'7;so[R] — £'^;So[R]| of the 
modes j,(3. This approximation will be satisfied if the ratio Jh/Jp ^ 1, as the 
separation of the modes is at least as large as Jp, see figure for an illustration. 
However, this approximation can have a much larger range of validity in practice, if 
the energy difference |i?7;So[R] — -E'/3:So[R]| is greater than or equal to the binding 
energy of a single ALE. For example, this approximation is clearly also valid provided 
Jh «C U, c.f. figure [3b. 

The diagonal terms a = /S in the second rhs.-term of (fTB]) give an effective 
renormalisation of the hopping rate Jh that appears in the kinetic energy term in 
the first line of the same equation. This assumes that the coefficients given by the 
trace over r are approximately equal for all R, which we have confirmed by evaluation 
of the corresponding terms from exact calculations in ID. We find maximum variation 
of the coefficients of the order of 15% even when two holes are very close (within 
2 sites), whilst for a more dilute gas, this approximation becomes even better. 
The renormalised kinetic energy can be interpreted as the kinetic energy term of the 
composite exciton, and in the tightly bound limit, \ Jp\, \ Jh\ U, reduces to the values 
found in sec. 12.21 

The solution offTSlis simplified in our case because the particles are non-interacting 
fermions, and thus the problem reduces to finding the solutions for a particle on a 
lattice with N impurities. This is detailed in [Appendix B[ We assume throughout 
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Figure 3. Diagonal and nondiagonal contributions of the last term of III6I I 
for Af = 3 and Jh/Jp = 0.1 for a 60-site lattice. (a) Jbo plotted as a 
function of U /\Jp\, computed from the diagonal elements. As ALEs become more 
tightly bound, the result approaches that obtained from perturbation theory, 
Jsff = "iJpJh/U. (b) Maximal effective coupling between the lowest channel 
7 and the first three excited channels /3i, /32, arising from the off-diagonal 
elements, max_R,„{|Jp E,0;(2:;^)?i/3fc ^ ± e")/(£;^.so [5 ± e"] " -B^;SO [«])!}■ 
Even for low U /\Jp\, nonadiabatic coupling will always add just small perturbative 
corrections, provided \Jh/Jp\ is sufficiently small. 



that the N impurities are sufficiently separated from each other that TV bound states 
exist. The critical separation at which bound states occur is dependent upon U /\Jp\ 
and the number of particles. For N — 2 bound states exist even when holes are on 
neighbouring sites provided U/\Jp\ > 2. 

For N — 2, The single-particle bound states are then the symmetric and 
antisymmetric superpositions p± (jB.3|) of two exponentially decaying bound states, 
centred at Ri and R2. These take the same form as the single ALE wavefunction, ([S]), 
but now with different energy E that depends on the separation R — |Ri — R2I. 
The lowest-energy two-atom wavefunction 0(ri, r2; Ri, R2) is then obtained from the 
Slater-determinant of p±. In principle, the choice of overall phase of these solutions can 
be made differently for different Ri , R2 . Without loss of generality we can choose each 
(?!)(ri, r2; Ri, R2) to be real, leaving us with a choice as to whether (/)(ri, r2; Ri, R2) 
should be symmetric or antisymmetric under exchange of hole coordinates Ri and 
R2. However, in order to minimise the non-adiabatic terms in (flB)) . we are required 
to choose these to be antisymmetric under exchange of Ri and R2, in addition to 
the requirement of antisymmetry under exchange of particle coordinates ri and r2. 
It then follows that in order to ensure the correct antisymmetry for exchange of 
hole coordinates Ri in the total wavefunction ^/'(r, R, i), [(|14p]. that C-y(Ri,R2,t) 
must be symmetric under exchange of coordinates R;. In the limit where the BO 
approximation is valid, we can then interpret the coordinates R^ as approximate 
centre of mass coordinates for the ALEs, which behave under this approximation as 
composite bosons. 

The lowest values of the potential E^.bo[Ri, R2] in ID are shown for various 
values of J7/ 1 Jp I in figurelH where the zero of energy is chosen so that Ea-BO [^1,-^2] 
as i? = |i?i — i?2| — > 00. For large U /\Jp\, interactions are very short-ranged, reducing 
to the values found in section [^21 As U/\Jp\ is decreased, the increase in mobility 
of the conduction band atoms leads first to increased interaction strength, and then 
as U/\Jp\ is further decreased, to longer-range interactions. These are shown in the 
figure around the parameter values where they become significant over 2-3 sites, which 
is approximately between U/\Jp\ = 2 and U/\Jp\ = 3. For finite Jh these interactions 
will be important as long as the energy of the COM-motion is not much larger than 
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the interaction energy. 




Figure 4. Born-Oppenhoimor curves for the interaction between two ALEs as 
a function of separation _R, for different values of U/\Jp\. For U/\Jp\ ~ 2 — 3, 
longer-range interactions are significant. For tightly bound excitons however, the 
Born-Oppenheimer approximation recovers the perturbative limit, in which only 
nearest-neighbour interactions are appreciable. The zero of energy is chosen so 
that all curves are to zero in the limit R/a ^ 1. 

For N > 2, we can always decompose the potential Ea-Bo[^ into a sum over 
contributions from different numbers of particles. 

(2) 

For example, for = 3 we can take the two-body potential E^.^q [Ri, R2] found 
for two particles (when the third is an infinite distance from the first two), and define 
the 3-body interaction [Ri, R2, R3] such that 

i<j 

By solving (fT5|) for = 2 and = 3 we can thus assess the importance of genuine 
three body interactions. Typical results for the ID case, in which we are primarily 
interested, are ploted in figure[5l They show two typical solutions for E^^^q [0, R2, R3] , 
and we see that only for low U/\Jp\ there are any appreciable three-body interactions, 
even at small distances. The same will be true for larger N provided that the density 
remains sufficiently small. Thus, for cold ALEs at sufficiently small densities, the 
BO-potentials in ^B)) can be decomposed into a sum of two-hole interactions, 

i<j 

There are two cases where the symmetry properties of the two particle 
wavefunction also transfer to the A^-body case. The first is the limit of two-body 
collisions, where the only significant interactions are those between each particle and 
its nearest neighbour, 

Ec^-Boim « ^I]max4^^o[R„R,]. (19) 

t 

In this limit, we can construct 0Q(r;R;) from a Slater determinant of single particle 
bound states at single holes or combinations of two holes, and the arguments presented 
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Figure 5. Three-body interactions between ALEs E^^/'^qIO, R2, R3] computed 
in the Born-Oppenheimer approximation curves for (a) U/\Jp\ = 2.1 and (b) 
U/\Jp\ = 3. These results should be compared to the two-body interactions 
plotted in figure |4] Even when the three ALEs are close together, the three-body 
interactions are dominated by two-body interactions, with three-body interactions 
becoming extremely small for U/\Jp\ > 3. This is still the case near the limit 
U/\Jp\ ~ 2, at which limit one bound state disappears for the case of ALEs 
existing on neighbouring lattice sites. 



for the two-particle case will generalise. The other case is that in which the holes 
are evenly spaced. There, the correct wavefunction is again the Slater determinant 
of single particle bound states centred at each hole, irrespective of the separation 
provided that the N bound states exist. This problem can be seen analogously to 
that of a periodic potential in free space, with the single particle bound states at each 
hole playing the role of Wannier functions for this periodic structure. In each of these 
cases, the optimal choice for the Born-Oppenheimer approximation is that (/'Q(r;R) is 
antisymmetric under exchange of hole or particle coordinates, and it follows that the 
wavefunction C^(i?i, . . . , Rn) will be symmetric under exchange of coordinates Ri, as 
discussed above for the N = 2 case. 

In either of these regimes, we then can write an effective Hamiltonian for the 
ALEs as composite bosons with creation operators b^, which in ID takes the form 

Hhc = -Jbo J2 ^l^i + J2 ^i^^^^+u Vi = ^i'Ll^o, Ri], (20) 

hi — b\bj, and Jbo is the renormalised effective hopping given by the first and last 
terms in (llOp . It can be approximately calculated from the solution to the single 
impurity problem {N = 1), (j)^{ri — i?i), as 

Jbo -JhY, <^7(^i - Ri)4>i{ri - i?i + 1). (21) 

ri 

In figure [5^, Jbo is plotted. 

In this section we have effectively shown how imbalance in the hopping rates 
can enhance the strength and range of interactions between ALEs. This will lead 
to the stabilization of diagonal order, i.e. crystaline order away from half-filling, 
which is discussed in detail in Sec. 13.11 Note that this limit is extremely relevant for 
ALEs, where Jp will typically be an order of magnitude larger than J^, and where the 
lifetime of an ALE is sufficiently long for longer range interactions to form a crystaline 
structure. 

2-4-. Interaction of ALEs: Numerical results 

Beyond the tight-binding and Born-Oppenheimer cases, we have also performed small- 
scale exact diagonalizations for two interacting excitons in ID, using parameters 
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ranging from U/\J\ = 1 to U /\J\ — 25 on up to 22 lattice sites. Strong imbalances 
in the hopping rates of holes and particles were considered as well, ranging from 
Jh/Jp = 1 to Jh/Jp = 1/40. All these calculations show only effective repulsion 
between the two excitons and a complete absence of bound states of four particles. 
A typical result is depicted in figure [6l where the scattering continuum and the 
dispersion relation of the anti-bound state resulting from weak repulsion can be clearly 
distinguished. 

In all cases in ID that we have analyzed, both analytically and numerically, 
we only find repulsive interactions, and thus no bound states. With increasing 
spatial dimensionality, it becomes generally more unlikely for bound states to form in 
scattering problems, and thus we also expect that this result will hold for arbitrary 
values of \Jh\/U, \Jp\/U in 2D and 3D. 
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Figure 6. Discrete energy eigenvalues for two interacting excitons on 22 lattice 
sites, plotted as a function of quasimomentum. Black crosses correspond to 
scattering states. Orange crosses form the dispersion relation of a repulsively 
bound state resulting from weak exciton-exciton repulsion. 



3. Many ALEs at low temperatures - exciton condensates and crystals 

Much fundamental interest in semiconductor excitons has stemmed from their 
predicted ability to undergo condensation in 3D or to form quasicondensates in 
ID. Due to their isolation from the environment, ALEs present an opportunity to 
investigate exciton condensation in a clean and highly controllable environment. In 
addition, the long lifetimes of ALEs, combined with long range interactions that are 
found in the case of unbalanced hopping give rise to other potentially interesting 
phases, especially an ALE crystal. 
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3.1. The ALE condensate 

Due to the effective bosonic nature of ALEs, we expect them to form a condensate 
characterised by off-diagonal Long-Range Order (ODLRO) at sufficiently low 
temperatures, in analogy to semiconductor excitons [2]. This should be true 
irrespective of whether the excitons are strongly or weakly bound. 



3.1.1. Continuum approach The mean-field approximation normally made for 
semiconductor excitons is a continuum model, and can be generalised to the lattice 
in order to obtain a description of the ALE groundstate at zero-temperature. In the 
continuum limit we write the Hamiltonian as 



H = 



E 

k 

1 



E 

klq 



qC|-qClCk 



-V^d^ 

2 "^q^k 



Vc 



q4+q4-q^lCk 



(22) 



where i?^ and are the dispersion relations for particles and holes respectively, and 
' ^q* ^'^'^ ^ denote the effective interaction potentials for particle-particle, hole- 
hole, and particle-hole interactions, respectively. The exciton creation operator can 
then be written as 



it|0) = 5]A„(k,k')44'|0), 

k,k 



(23) 



where the index a will specify the quantum numbers for the exciton. The key 
observation is that the single-exciton state of lowest energy will have the lowest 
possible bound-state level n = 0, and COM-momentum K = w. Here, w denotes the 
difference in position between the maximum of the valence band and the minimum of 
the conduction band. It can either be or half a reciprocal lattice vector. If w = 
(as was the case for the ID ALEs in the previous section) the system has a direct hand 
gap. Whereas for w ^ it has an indirect hand gap (the 2D and 3D realization of 
ALEs are of this form) . 

In the following, the dispersion relations of the valencce band, is shifted by — w. 
Then the single exciton operator of lowest energy with a = (w, 0) has the form 
^oo|0) — X^k ^(^)'^k^-klO) with ^(k) — Ao,o(k, — k). Such a state corresponds to 
an effective pairing of electron at k and hole at w — k (when the notational shift of 
the hole-quasimomenta is reversed) . The equivalent of this pairing in an approximate 
many-body groundstate emerges from an ansatz first proposed by Keldysh and Kozlov 
[T]. There, the ground state of the system is assumed to be a coherent state for this 
single-exciton operator A\ g , 



V\Q) = exp 



^o,o) 



\{{uv. + v^c\dl^)\Q), (24) 



;^^(k)) wk 



(25) 



The resultant state has the same structure as the BCS ground state of weakly 
correlated electrons in a superconductor. At finite density the coefficients A{k) 
in ([M]) and (P5|) arc no longer the single-exciton ground state wavefunction. Instead, 
the coefficients Uk and Wk have to be determined in a self-consistent manner 
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(see I Appendix C\ for background). In mean-field theory and using the constraint 
of fixed density these are found to be 

^k=-E^'-'^^' «^^ = ^Efi-|tl (26) 



where 



a = E^ + E'i + 2Y,^^-^^^l (27) 
1 

Sk = V, + \{Vi + V>^-V^^V^) (28) 

A^ = -E^i-k"Fl (29) 
1 



+ (30) 

In this regime the Hamiltonian ([221) can be diagonahsed with new fermionic 
quasiparticle and quasihole operators Cj^ and D^^ [see (|C.1I) ] to yield 

Hmf = J2 [eL'^cIc^ + , (31) 

k 

with the expressions for the new quasiparticle dispersion relations E^^ given 

in [Appendix C| ([CJl) . ([CJll . 

As a consequence of this ansatz for the groundstate, the lowest energy solution 
(k, n) = (0, 0) of the single exciton can be recovered from the above mean-field 
equations in the limit of vanishing density Uf^^ —^ [41| . 

The condensate thus exhibits off-diagonal long-range order, with the associated 
dissipationless transport properties. For excitons, however, these involve only 
dissipation-free transport of excitation energy and momentum, not of mass or charge 
transfer. 



3.1.2. ALE condensate on a lattice We can adapt Hamiltonian (I22|) - and the mean- 
field results derived from it - directly to the case of ALEs, where k now denotes 
lattice quasimomenta in the first Brillouin-zone (B.Z.) instead of momenta in the 
continuum case, and we set Vq = 0, Vq = 0, Vq = —U. Using Hamiltonian and 
applying (j26p - (j30p . the condensate at zero temperature can be approximated by a 
BCS-type groundstate 

\^)=X{{n^ + v^cldl^m (32) 

k 

where Uk and Wk must obey the gap-equation and the density constraint [HI [42l [43] : 
U I 77^^, ne^ = —l ( 1 _ _ ) (33) 

where 



a = - E('^P + cos(fcrfa) - flex + 2Unex, (34) 



E^^Ja + {2APf, (35) 
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and Jp denotes the hopping rates of conduction band atom in the d-th direction. 
As the conduction band is assumed anisotropic (see Sec. [2]), ^ Jp,Jp wiU hold 
generally, with < 0, Jp,Jp > [57]. Consequently, the conduction band has its 
minima around (±7r/a, 0, 0), and thus w = (0, tt/o, ir/a). Contrary to the situation in 
ID, the system thus has an indirect gap, and at zero temperature the ALEs condense 
into a motional state with nonzero quasimomentum. 

Evaluating (|C.6P and (|C.7|1 for the given model of ALEs yields the dispersion- 
relations of the new quasiparticles and quasiholes respectively: 

= liK- + E^) , = 1 {E^^ + E^) (36) 

where 

D D 

EP = -2Y,JpC0s{kda), Et^~2Y,JhCos{kda). (37) 

d=l d=l 

For two-band implementation of ALEs considered here, where ^ Jh in 2D and 3D, 

the quasiparticle dispersion relations are different, E^"^ ^ Thus they exhibit 

minima at different points in the Brillouin zone. 

As shown in [JTJ this BCS-formalism can interpolate between the weak- 
coupling U \Jp\, \Jh\ and strong-coupling U ^ \Jp\, \Jh\ regime, provided the 
Hartree-Fock corrections to the chemical potential are included. For strong coupling it 
also becomes exact in the low-density regime (n^x ~ 0) and provides a good qualitative 
description of the groundstate for increasing density up to the case of maximal filling. 

In Sec. m we will discuss measurements that can be made on exciton condensates, 
including determination of the condensate fraction and pairing correlations. 

3.2. The ALE crystal in ID 

The long range repulsive interactions that we found for large hopping imbalances 
Jp ^ Jh in Sec. 12.21 suggest the possibility of forming an ALE crystal for densities 
away from half filling, n^^x 7^ \- This situation is directly related to recent research on 
Fermi Hubbard models with imbalanced hopping rates for the two species [40i i44) . 

The results of the Born-Oppenheimer approximation in Sec. 12.31 suggest that the 
effective interactions have maximal range and strength for loosely bound ALEs at 
U/\Jp\ ~ 2 - 3 for large hopping imbalances. For these parameters, naive inspection 
of figure [4] suggests that interaction effects will favour Diagonal Long Rangle order 
(DLRO), i.e., a crystal, and suppress pairing maximally from n^x ~ ^ upward, 
provided that the kinetic energy of the ALEs is sufficiently small. We investigate the 
possibility of such a crystal using imaginary time evolution on matrix product states 
(with a timc-dcpndcnt DMRG algorithm) [45l |46l [47l |48] to compute the groundstate 
Hamiltonian ([3]) at different densities. From the numerically calculated groundstates, 
we obtained density-density- and pairing correlations for four different densities, with 
nex = M if' i§ ^^'^ if' ^ 96-site lattice for different values of U/\Jp\ and Jh/ Jp- 

In bosonised theories [49 , the Luttinger parameter K determines the algebraic 
decay of correlation functions in the thermodynamic limit at large distances (true 
long-range order of any type being absent in ID for continuously- valued fields). We 
expect the pairing correlation function {cl.dl.dyCy) to behave as Ap/\x — y\^'' , and the 
density-density-correlation function {{uxny) — {nx){ny)) to behave as Add cos(fc/(a; — 



Atomic lattice excitons: from condensates to crystals 



15 



y))/\x — y\^o" for \x — y\'^ 1, where kf = nuex denotes the Fermivector. Ap, Add, 
Kp and Kp)p) were determined from fitting to these functions. 

Bosonisation predicts Kp/p, = 1/Kp [JS]. The vahdity of bosonization will 
shift to decreasing wavelengths, however, as Jh/Jp decreases. This results from the 
increasingly different slopes of the dispersion-relation linearisation around the Fermi- 
points. Excitations that are still weak - i.e. long-wavelength - for the more mobile 
conduction-band atoms, will be in the non-linear regime for the slower atomic holes. 
The validity of bosonization for any finite-size system is thus limited by the length- 
scale set by the less mobile holes, and we expect our simulation results to exhibit 
noticeable deviations from bosonization predictions below some value of Jh/Jp- 

While we averaged out fc/-oscillations superimposed on the algebraically decaying 
pairing correlations (c.f. [50]) to obtain reliable fits for Kp, this is more difficult for 
Kpp in the density-density-correlations. As a result, we obtain large deviations from 
the bosonization prediction Kpp = 1/Kp even for Jp = Jh, where bosonization 
is most reliable. Despite this difficulty, we can identify clearly different regimes of 
behaviour, dependent on the hopping imbalance and ALE density. For illustration, 
the powers of algebraic decay obtained from the fit are depicted in figure [7] for the 
pairing correlations, while figure [5] shows representative examples of density-density 
correlation functions in different parameter regimes. Generally we observe in our 
results that superfluid order is dominant for equal or moderately imbalanced hopping 
rates, but becomes more and more suppressed as ALE density and hopping imbalance 
increase. Complementary, density-density-correlations start decaying slower and 
become the dominant order. As is to be expected in ID, this superfluid-to-crystal 
transition is continuous. Crystaline behaviour is especially pronounced in the region 
predicted from Born-Oppenheimer calculations, > i and U/\Jp\ = 2 

— 3, where density-density-correlations decay only moderately (figure [8]) and ALE- 
interactions suppress pairing very strongly (figure [7]). 

For higher dimensions, the subject of crystal formation close to half-filling for 
general values of U/\Jp\,U/\Jh\ has recently been addressed using both numerical 
DMFT and analytical mean-field analysis 01]. The results are directly applicable 
to loosely bound ALEs. Both DMFT and standard mean-field calculations show 
the ALEs to form a superfluid condensate for moderate hopping imbalances, and to 
phase separate into checkerboard crystal and superfluid below some critical, density- 
dependent value of Jh/Jp (c.f. figure 1 in [44j). For a parameter regime analogous 
to the one considered here (large hopping imbalance, loosely bound ALEs) analytical 
mean-fleld calculations additionally predict the possibility of pure crystal phase at 
noncommensurate densities (c.f. flgure 2 in [44]). As the DMFT calculations do not 
yield such a phase, the authors question the validity of this result. Whether ALE- 
interactions can thus suppress superfluidity at large hopping imbalances sufficiently for 
the formation of a stable noncommensurate crystal phase in higher spatial dimensions 
is not fully resolved yet. 

As a final comment, note that the Hamiltonian ([3]) with Jp ^ Jh used here 
to describe ALEs, has been studied previously in ID [40] using Luttinger liquid 
theory. There, the authors showed that increasing hopping imbalance can lower 
the value Kc, the exponent controlling the algebraic decay of correlation functions. 
For |x — y| ^ 1, density-density (DD) correlations decay as ((rix^-y) — {nx){ny)) ^ 
cos(kf{x — y)))/\x — y\^'', with Ux = cl.Cx + d\.dx, whereas pairing-correlations decay as 
{cldldyCy) ~ l/\x — y]^^^" [49]. As Jh/Jp decreases, Kc can decrease below 1. Then 
crystaline ordering becomes dominant over superfluid ordering, which is dominant for 
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Figure 7. Power of algebraic decay of pairing correlations plotted as a function 
of U/\Jp\ for Jp/Jh = 1 and Jp/Jh = 0.025, and for four different densities 
on a 96-site lattice, (a) For Jh/Jp = 1, pairing correlations decay moderately 
faster with increasing density and decreasing U/\Jp\. (b) For strong hopping 
imbalance, pairing correlations decay much faster with decreasing U/\Jp\ as ALE- 
ALE interaction increases in range and magnitude, which become more significant 
as Uex increases. The decay is especially strong for n^x > 0.3 in the regime 
U/\Jp\ 2 — 3. From the Born-Oppenheimer theory this result is expected, 
see e.g. figure l4l as for these parameters and Jh/Jp = 0.025 the next-nearest- 
neighbour interactions between ALEs is still dominant over the kinetic energy 
of the COM-motion of the ALEs. As a result, for nex = 0.31, close to the 
ideal crystal filling factor Uex = ^, the decay also shows the steepest increase 
as U/\Jp\ decreases. The enhanced next- nearest-neighbour interactions in this 
region strongly suppress pairing order, and thus for the ALEs crystaline order 
dominates in this regime. 
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Figure 8. Representative examples of real-space density-density correlations for 
Jp/Jh = 1 (solid lines) and Jp/Jh = 0.025 (dotted lines) plotted from the center 
of a 96-site lattice. Here, 1549,1/ = (nigriy) — {n4Q){ny) is shown at four different 
ALE -densities for U/\Jp\ = 3. At all densities, both amplitude and rate of decay 
of these DD-correlations are significantly improved by hopping imbalance over the 
case of Jp = Jh- The overall magnitude of the correlations increases with Uex- 
Meanwhile, pairing correlations start to decay rapidly in this regime of U/\Jp\, 
see figure[7] Taken together, this identifies the regime of \ Jh/Jp\ <S 1, U/\Jp\ = 2 
— 3 and riex > 0.3 as the region where the crystaline characteristics of the ALEs 
are most pronounced. 
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Kc > 1. The value of Jh/ Jp where = 1 wih generally depend on the exciton density 
Uex- Based on the value of a phase diagram can be drawn, c.f. figure 1 in [40] . 

3.3. Phases of tightly bound ALEs 

As ALEs map onto an attractive Hubbard model with imbalanced hopping of particles 
and holes, there are also some previously known results for case of isotropic hopping 
Jp = Jp ~ Jp = Jpj yielding the phase diagram in the tightly bound limit 
\Jp\/U, \Jh\/U ^ 1. Here ALEs can be described as hard-core bosons with hopping 
rate Jeff and repulsive nearest-neighbour interactions Veff. For a regular square or 
cubic lattice, mean- field calculations (c.f. e.g. [35]) predict either a checkerboard 
'crystaline' Charge-Density- Wave (CDW) or a superfluid order at maximal density 
nex = 5, depending on the value of 14/// Je//. The CDW is characterised by (DLRO), 
as was the exciton crystal studied in the previous section. Away from maximal 
filling, < 5, mean- field calculations predict a supersolid phase - the coexistence 
of superfluid and crystaline order - for a broad range of n^x (c.f. figure 4 in [39 ). 
Beyond this range of densities, pure superfluidity is predicted. However, Quantum 
Monte-Carlo simulations [STl |52l [Ml [54j show the prediction of a supersolid away 
from commensurability to be inaccurate for regular lattices. Depending on Vef f / Jef f 
and n^xi tightly bound ALEs are either superfluid or undergo phase separation into 
spatially disjoint subsystems, one exhibiting superfluidity, the other checkerboard 
crystaline order [^. Both Dynamical Mean Field Theory (DMFT) and standard mean- 
field calculations suggest ^44j that when Jp = Jh (i.e. Veff /Jeff = 1) tightly 
bound ALEs will be superfluid at any density. One way to increase the value of 
^eff/Jeff — {Jp + JD/C^JpJh) and thereby attain phase separation beyond some 
critical n^x is to have Jp ^ Jh (c.f. [H]). The realization of ALEs proposed here will 
always have a signiflcant imbalance, as the conduction band hopping rate typically 
surpasses that of valence band by an order of magnitude. Another possibility to 
increase the ratio Vef f / Jef f is to introduce additional nearest-neighbour interactions 
in the Hamiltonian ([3|) (see [39]), e.g. by considering fermionic dipolar molecules 
instead of atoms. 

In summary, while for sufficiently low Uex, tightly bound excitons always 
exhibit just superfiuid order, exciton-exciton interactions can enhance density-density 
correlations. For larger values of n^x and with hopping imbalance between particles 
and holes, they can result in phase-separation of the system into a checkerboard ALE- 
crystal and an ALE superfluid. 

4. Probing atomic lattice excitons 

Optical Lattice experiments provide a range of measurement possibilities that we 
can use to determine the characteristic properties of the many-body states that are 
produced. Three important examples of this for probing ALEs are Lattice modulation 
spectroscopy, RF Spectroscopy, and noise correlation measurements. 

§ On triangular - and thus frustrated - lattices, Monte-Carlo simulations |56| have shown that a 
supersolid phase can appear however. 
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Jp\+\Jh\ 



Figure 9. (a) Sketch of ALE phase diagram as a function of hopping imbalance 
and interaction strength relative to kinetic energy(after figure 2 in 1441 ). 
(b) Sketch of phase diagram for tightly bound ALEs at fixed density. SF: 
Superfluid, PS; Phase separation, CDW: checkerboard Charge-Density Wave, 
ncCDW: Charge-Density Wave at noncommensurate density, (a) Mean-field phase 
diagram for ALEs at fixed density close to riex = ^ ■ When the hopping imbalance 
Jp I — [ Jj, I / 1 Jp I -I- 1 Jji I is large enough, the superfluid phase-separates into superfluid 
and checkerboard-CDW. The separating line shifts with density. A pure non- 
checkerboard CDW-phase might also be possible in a small parameter range. For 
large U /{\Jp\ + \Jh\), ALEs become hard-core bosons with ranged interactions, (b) 
Sketch of a qualitative phase diagram for tightly bound ALEs. Doped away from 
the checkerboard CDW phase at n^x = \, ALEs phase-separate into superfluid 
and checkerboard CDW, if V^f f / J^f f is largo enough. This ratio can be tuned 
by changing / Jp or by introducing longer-ranged interactions between particles 
and holes. 



4-.1. Exciton detection via lattice modulation 

Periodic modulation of the lattice depth has become a standard way to probe the 
excitation spectrum of many-body systems on an optical lattice [55]. It is usually 
applied to systems in the groundstate. If the frequency of the modulation matches 
that of an excitation of the system, it drives transitions to higher-energy states. Such 
a transition is then detected by ramping down the confining potentials/interactions 
adiabatically, thereby transferring all potential energy into kinetic energy of the 
expanding atom cloud, which is then measured. This method has also been used 
to drive transition from excited states to lower energy states recently, in order to 
determine the energy of pairs of atoms bound together by on-site repulsion |17j . 

Qualitatively, probing for excitons via lattice modulations should show a 
pronounced decrease in the systems total energy around a frequency given by the 
bandgap minus the binding energy of the exciton. 

4-2. Measuring the exciton condensate fraction via RF- spectroscopy 

In order to demonstrate particle-hole pairing and in particular to detect the 
macroscopic occupation of the exciton groundstate, it is possible to employ a range of 
measuring techniques that have been developed for cold atoms experiments. At zero 
temperature, with no thermal excitation of collective exciton modes with K 7^ w, the 
pairing amplitude is directly proportional to the condensate fraction. It can be 
probed by the application of RF-spectroscopy, following the example of [17] (c.f. [5H] 
for an experimental application). An RF-pulse with frequency lorf is used to couple 
the conduction band atoms, at energy ujc, to another internal state with energy W3, 
with a detuning 5 — ujrf — (^3 — Wc), c.f. figure fTOl Choosing (5 > 0, the particle- hole 
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pairs are broken up in this process. It is assumed that this state is still lattice-confined, 
with single particle energy /^3,e//j where e^. is the dispersion relation of 

the state on the lattice and fJ-s.eff is the effective chemical potential, i.e. the chemical 
potential including any mean-field shift from interactions. 

The RF-pulse with frequency lorf has an amplitude f2 that can be assumed to 
be slowly varying on length scales of the lattice. The transfer rate into the third 
state as a function of the RF-detuning, I{S), can be calculated by generalizing the 
results in [57^ to the case of ALEs with anisotropic conduction band hopping (and 



therefore E^ ' ^ E^ '). The result (at zero temperature) then becomes 



I{S) = -2n\n\' ^ vlSi^i'^ + A + 4^))0(^, 



(38) 



Here, A = /i3 — /i^ — <5, and S ~ liJrf ~ (^3 ~ "-^c) denotes the detuning of the RF-pulse 
from the transition frequency between the two internal states, /ic is the chemical 
potential of the conduction band atoms, for which relation /ic + M/i = Mex holds. This 
expression is then valid generally for any Jp and J/i, assuming positive detuning. 

The spectrum described by (I38p has a gap, given by the minimum value of the 



(5- function argument. Its value is 8g 



J3) 



' mink E^ — jic, assuming the additional 



internal state is initially unpopulated and does not scatter off the other states, i.e 
Hs.eff = f^s — 0- If nex is known, A^ can be determined selfconsistently from 6gap 
using ((35)) . 



UJr: 



a) 



third level 

- conduction band 

- valence band 




Figure 10. RF spectroscopy for an ALE condensate, (a) Excitation with RF- 
coupling Q and detuning 5 between conduction band atoms and a third internal 
state, (b) Transfer rate |/((5)/27rf!^ | over detuning 5 for U/\Jp\ = —8, U/\Jh\ 
-40, rie 

The kinks result from van Hove singularities. 



: 0.3. The gap in the spectrum Sgap is given by c'q ^ + miuj^ E^"' — /x^ 



4-3. Detecting crystal structure and pairing correlations via noise spectroscopy 

We can detect both crystal structure and pairing correlation of ALEs using atom 
shot-noise measurements, a technique that was proposed in [35] and has been 
demonstrated in the laboratory (3^1 ISZ] ■ The former can be identified through second- 
order correlation functions of valence-band atoms, and the latter from second-order 
correlation functions of atoms in two bands. Each of these can be obtained from 
fiuctuations in the density profile of the atomic gas. 

For the ALE crystal, we expect holes to be anticorrelated with each other around 
the Fermi-edge. As the operators bt for valence-band atoms in are related to the 
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hole operators through 6^ ~ "^-k' thus hkewise expect an anticorrelation peak 
around the Fermi-edge for the valence-band atoms. In the case of the ALE-condensate, 
which is characterised by the pairing of particles and holes at opposite quasimomenta, 
pronouced anticorrelation of conduction- and valence-band atoms should be visible at 
equal quasimomentum around the Fermi-edge, i.e. " antipairing" . 

To obtain the density profiles experimentally, and from them the fluctuations, 
the Brillouin zones need to be imaged, i.e. lattice quasimomentum k needs to be 
mapped to real space position R(k) on the detector. This is achieved by ramping down 
the lattice sufficiently slowly to keep the atoms within their respective bands whilst 
preserving their quasimomentum (c.f. [13j). Valence-band atoms then occupy the first, 
and conduction-band atoms the second Brillouin zone (c.f. figure [Tib and [T5b) . If 
density fluctuations are just limited by shot noise, density fluctuations of conduction 
and valence band atoms will be correlated according to the connected correlation 
function 

t;,„(R,R') = (n^(k)nV)) - ('^R(k))(^R(k'))' (39) 
whereas the density- fluctuations of valence-band atoms are correlated amongst 
themselves according to 

t;.„(R,R') = (n^(k)nV)> - (^^R(k)>("R(k'))- (40) 

Here n^^^) = 4(k)CR(k) and n^^^,) = 6R(k')^R(k')- 

For the case of ALEs in ID we used the numerical algorithm employed in 
section 13.21 to compute both correlation functions in different parameter regimes. 
These are plotted in Figs. [11] and [T^l which when contrasted display the crystal 
and superfiuid characteristics of the system in different parameter regimes, figure [TTk 
shows the momentum-space density profile for excitons in the first and second Brillouin 
zones for the example of a parameter regime where we expect a crystal to form. The 
fluctuations that would occur in the experimental measurement of this density proflle 
can be used to compute the second order correlation functions ([5^1) and (|T0|) [55] . 
Figure [TTb depicts Gcv{R,R'), that results from Brillouin-zone resolved measurement, 
for different system parameters. It shows how antipairing is suppressed when nex 
increases and U/\Jp\ decreases. Figure [T^ shows the periodic dips at G ± UexT^ja 
in GcviR,R') (where G is any reciprocal lattice vector). This can be measured 
from the noise correlations measured from real momentum distributions, as is done 
in switching off the lattice suddenly, analogous to previous experiments [SHI [57] . 
The periodic anti-correlations indicate the crystal, and their locations relative to the 
large positive correlations reflect the reciprocal lattice vector of the exciton crystal. 
Dips appear here because we are computing noise correlations for fermions |35j . 
Taken together, the figures further illustrate the analysis of section [3721 They show 
that pairing and density-density correlations are complementary, with pairing being 
strongly suppressed and density-density correlations around the Fermi-edge enhanced 
in the limit of strong hopping imbalance and n^x ~ U/\Jp\ = 2 — 3. This 
indicates onset of the exciton crystal. Paring increases and the dips in the density- 
density correlations vanish when Uex decreases or Jh/Jp increases, corresponding to 
the dominance of superfiuid character in the system. 

For higher dimensional systems, the existence of a pure crystal phase is not 
resolved, though crystal and superfiuid order can coexist (see sec. 13. 3|) . When only 
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Figure 11. (a) Density profile for quasimomentum states in the first two Brillouin 
zones for ID ALEs, calculated for a 40-site lattice in the crystal regime (riex = ^ , 
U/\Jp\ — Jh/Jp = 0.025). (b) Pairing anticorrelations at fc = fc' between 
conduction- and valence-band atoms in different parameter regimes on a 40-site 
lattice. With parameters such that we expect to observe an ALE crystal (solid 
line, n^x = U/\Jp\ = 4, Jf^/Jp = 0.025), antipairing is strongly supressed 
as compared to a system outside this regime (other lines), with a minimum at 
the Fermi-edge, Uex't^/o. away from the minimum of the Bloch band at ir/a. For 
equal hopping (dotted line, n^x = ^g, U/\Jp\ = 4, Jh/Jp = 1) antipairing is 
strong and dominant (see figure [121 . At lower densities (dashed line, (n^x = 
U/\Jp\ =4, Jh/Jp = 1) antipairing is still stronger than for the cystal, with the 
minimum shifted due to lower density. For increased attraction (dash-dotted line, 
riex = -jj, U/\Jp\ = 10, Jh/Jp = 1) antipairing is slightly increased over the 
crystal 
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Figure 12. Second order correlation function for valence band atoms, computed 
on a 40-site lattice. These represent results that can be obtained from noise 
correlation measurements, in which the crystal structure of ALEs appears. In 
the appropriate regime (black line, Uex = U/\Jp\ = 4, Jh/Jp = 0.025), 
dips appear at itnea;7r/a away from the location of reciprocal lattice vectors 
corresponding to the optical lattice. These signalize the formation of the ALE 
crystal, as holes become localized through atom-mediated long-range repulsions 
(see figure |4]| , and occur at momentum values corresponding to the reciprocal 
lattice vector of the exciton crystal. Such a result is representative of the crystal 
regime. If we increase U/\Jp\ to 10 (not shown), the dip visibility decreases only 
slightly, and the crystal structure still appears. Increasing Jh/Jh to 1 (grey line) 
causes the crystal and thus the dips to disappear, as does decreasing n^x to ^ 
(not shown) 
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superfluid order is present, the correlations between conduction- and valence-band 
atoms Q(.vi^^'^') can be expressed through the mean-field coefficients u and v. 

R(k) = R(k') 
R(k) ^ R(k') ■ 

An example plot of the column-integrated correlation function contained in the shot 
noise is shown in figure [T5b . Considering that the correlation at k = k' equals 1*^^^ 



""R,(k)'^R(k) 



(41) 




First 
Brilloui 
Zone 




R{5.o) 



Figure 13. (a) Calculated density profile of conduction and valence band atoms 
as they would appear after ideal time-of-flight imaging. Dark areas denote high 
atom density. The third dimension is column-integrated. The first and second 
Brillouin zone are mapped to real-space positions for the condensate ground state 
with n^^ = 0.3, U/\J^\ = -8, J2/Ji = J|/Ji = -0.5, U/\Jh\ = -40. (b) 
Calculated example of noise correlations between conduction and valence band 
atoms, at positions R-(k) = R(k'). Here R(k) is taken along the dotted line in 
figure a, and R.(k') is taken along the solid line. The anticorrelation of conduction 
and valence band atoms is most pronounced around the Fermi-edge of atoms and 
holes. As Jp < 0, Jp, Jp > 0, the Fermi-edge in figure a lies along the k2 axis. The 
anticorrelation appears reduced because of density-integration along R(0, 0, fcs). 



in the mean-field theory, this approach provides an alternative for determining the 
condensate fraction via the definition of A^. 



5. Exciton Formation in an Optical Lattice 

5.1. General remarks 

One method to prepare excitons in an optical lattice would be to create a band 
insulator of spin polarised fermions in the valence band, and then to excite atoms 
to the conduction band, using, e.g., a Raman transition. This would be most directly 
analogous to the method by which excitons are produced in the context of solid state 
systems. 

However, in the context of Atomic lattice excitions, it would be desirable both to 
have control over the number of excitons produced, and to have a means of producing 
these excitons in the lowest possible many-body energy state. In contrast to excitons 
in the solid state, the weak coupling of Atomic lattice excitions to dissipative processes 
means that excitons prepared with a significant energy above that of the groundstate 
will not be naturally cooled to this state. As a consequence, particle-hole pairs cannot 
be created in an arbitrary manner if the ground state is the desired outcome. 
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5.2. A preparation scheme 

One scheme that we have investigated in detail for its efficiency in preparing a 
many-exciton state close to the condensate groundstate is a melting scheme adapted 
from f59l. Through adiabatic tuning of the system parameters, a chosen filling factor 
of ALEs can be prepared at essentially T = (this is in the same spirit as other 
adiabatic preparation proposals, such as [60]). 

We begin with a spin-polarised Fermi gas loaded into an optical lattice with an 
applied superlattice, which gives certain lattice sites an offset energy V . The initial 
configuration of atoms is loaded so that the lowest motional state in each lattice site 
is occupied, as well as the first excited state in wells with lower energy, as shown in 
figure [Hb. This can be achieved, e.g., using the methods discussed in [551 [CTl W^. 
The final filling factor is chosen by the periodicity of the superlattice with respect 
to the basic lattice period. In each case, we will want to generate at some point an 
interaction energy U between atoms in the first excited motional level and atoms in 
the ground motional level. 

This can be achieved, e.g., by introducing an off-resonant Raman coupling 
between the initial internal state of the spin-polarised fermions and a second internal 
state. If the detuning is chosen so that atoms in the excited motional level in the 
initial internal state are coupled to the second internal state near to resonance with 
the first motional level, then the coupling will be much further detuned for the atoms 
in the initial internal state in the lowest band than for those in the upper band. This 
will produce a different admixture of the second internal state for atoms in each band, 
and the non-identical internal states will allow s-wave interactions between atoms in 
different bands. 

A simple scheme for the adiabatic formation of excitons is now the following: (i) 
Atoms are removed from the lowest motional level in sites containing two atoms, by 
transferring them to a different internal state and removing them from the lattice 
(see figure [T4b-I). This is possible because the resonant energy of the transition to 
a different internal state is shifted by differences in interaction energies, and because 
it is possible to address different bands using spin-dependent lattices for the initial 
and final states of the transfer [59]. (ii) The interaction U is "switched on", (iii) 
V is adiabatically decreased to zero, effectively melting the preformed excitons and 
forming the desired exciton gas (see figure iMh ) . 

Until the last step, the state is protected by a gap of order V, but the last step 
must be performed adiabatically with respect to the exciton tunnelling rate, which for 
Jp,Jh ^ U can be estimated as J ex ~ "^JpJh/U . 

These timescales can be improved upon by slightly modifying the above scheme. 
Instead of removing atoms initially in step (i), V is adiabatically decreased to zero, 
forming the ground state for delocalised fermions in the first Bloch band, as shown 
in figure [Hb-II. If the lattice depth is then suddenly increased to a large value where 
Jp, Jh ~ and atoms then removed as in (i), then we preform the ground state of hard- 
core excitons. If we then switch on the interaction U and adiabatically decrease the 
lattice depth, we obtain the desired exciton ground state. Note that in this scheme we 
must take care of the sign of the effective hopping for excitons, which in some cases 
requires an additional complication. The atoms that we delocalise in the first step 
should have the same hopping rate as the composite exciton. Thus, for Jp, Jh < 0, 
and U > 0, we obtain J^x > 0, and therefore atoms should first be delocalised in 
the lowest band, then transferred to the upper band whilst the lattice is very deep. 
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This requires atoms to be initially loaded in two internal states to allow for double 
occupation in the lowest band. 

The advantages of the latter scheme are clear from figure 1151 where we show the 
energy of the final state in ID as a function of ramp speed in the final step. We produce 
a final state with 15 excitons on 60 lattice sites, with U/\Jp\ = U/\Jh\ = 20, by ramping 
V and then J on a timescale given in units of Jg^^i ■ These results were computed using 
a time-dependent DMRG algorithm [45l |46l |47l |48] , and show clearly that ramping on 
a timescale much smaller than that implied by 2JpJh/U produce states with energies 
very near the ground state (indicated by the dashed line in the figure). The main curve 
shows the result for adiabatic ramping in the last step beginning from the grounds 
state of hard core atoms in single atoms, and the inset shows the energy of a single 
atom state after a similar ramp. 



a) 



b)® 
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Figure 14. Melting scheme to obtain an ALE condensate, (a) Initially the 
ground state is formed in the presence of a superlattice, giving site offsets V. 
Two alternative melting schemes, I and II, proceed from there. I: (b) ALEs are 
preformed locally by removal of the atom in the lower motional state on doubly 
occupied sites, c) Melting the state by ramping V down adiabatically yields 
delocalised ALEs. II) (b) Groundstate of delocalised fermions in the upper band 
is formed, melting the initial state by adiabatically decreasing V. (Depending 
on sign of the hopping for ALEs, this must be done in the lower band instead, 
see text for details) (c) The resulting state is frozen by ramping up the lattice 
suddenly, and atoms in the lowest band are removed from doubly occupied sites. 
This automatically pairs each atom in the upper band with a hole in the lower 
band, in the localised state expected in the limit U/J — > oo. Ramping the lattice 
down adiabatically then yields ALEs close to the ground state. 



6. Summary and Outlook 

Excitons are composite objects of fundamental interest in the context of semiconductor 
physics. However, the theoretical models presented for these systems relate even more 
closely to the situation we discuss in optical lattices, making ALEs an important 
testbed for study of the most interesting properties of these objects. In particular, 
the availability of techniques for state preparation and measurement in optical 
lattices provide tools to study excitons that strongly completment those available in 
semiconductor systems. ALEs can also be realised in parameter regimes that strongly 
contrast with those available in semiconductors. In this way we not only obtain the 
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Figure 15. Results of the Melting scheme II, as depicted in figure [T4l Main graph 
shows final state energy Ej/Jp obtained from numerical calculation of decreasing 
U I J for ALEs that were preformed utilizing the groundstate of delocalised atoms 
(see figure [Tib-II. c). Parameters are: 60 sites, 15 ALEs, V l\Jp\ = U/\.Jh\ = 20. 
Exact groundstate energy is indicated through a dashed red line. The inset shows 
Ef /J for the preparation of the groundstate of free atoms in the upper band from 
the intial band insulator (see figure [Mb . b-II), for identical parameters. Exact 
groundstate energy is indicated by dashed red line. 



possibility for condensation of these objects, but also for preparation of additional 
phases, including an exciton crystal. Optical Lattices also provide the possibility to 
directly reduce the dimensionality of the system, making the important case of ID 
excitons directly accessible. 

There are several open paths for the study of ALEs. Firstly, many of the results 
we have discussed in this work, particularly pertaining to the exciton crystal could 
be generalised to 2D and 3D. The structure of degenerate excited bands ^27] could 
also be introduced as an extra element in the higher dimensional context. In addition, 
the study of composite objects could be developed in the direction of multiple bound 
ALEs, in the case that, e.g., attractive interactions are generated between excitons. 
Other possibilities also arise to study systems analogous to excitonium (see [2], ch. 
10), a dynamically created collection of excitons that exists in the limit where the 
gap between valence of conduction-band is small compared with the atom-atom 
interactions. 

As a final comment, it is clear that there is a strongly connection between ALEs 
and superconductivity or superfluidity of fermions via the attractive Hubbard model, 
which we obtain here for interacting particles and holes. It would also be possible to 
consider the case of repulsive interactions between particles and holes, in which case 
one could obtain repulsively bound ALEs, in the sense that the two particle bound 
state appears above the scattering continuum. The model describing this system 
would be a Hubbard model with repulsive interations, which may lead to interesting 
analogies with d-wave pairing states. 
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Appendix A. Solution for the Single Exciton 

As for the COM-motion, ^ with [/ = corresponds to a Schrodinger equation for 
a free particle on a lattice. That equation is resolved by the Green's function with 
discrete translation invariance, 

r ( M dfc exp[ifc(r-/)] 

Go\r " r ) ~ a I —— — — r, — ^ ■ (A.l) 

^ ' J_^/^2t: E + 2Jcos{ka) + iT] ^ ' 

For U > - the single particle motion with an impurity at lattice-site r = - the 
Kroneker-delta in ([6]) allows an easy, explicit resummation of the recursive relation for 
the full Green's function: 

One can read off from this expression that the full Green's function has a single bound 
state pole for each value of K with energy Ek, where the denominator of the second 
term vanishes, i.e., for U^^ = Go(0). The numerator yields directly the unnormalised 
bound state wavefunction in relative coordinates, pr = CGo(r), which is given by 

_ dfc exp[ifcr] 

where C is the normalisation constant. This is depicted in figure[2]a,b. The dispersion 
relation E of the exciton, which forms the exciton solution in figure [2] is determined 
by the implicit equation 

-1 _ dfc 1 

^ 2^^ + 2Jcos(fca)' ^ ' 



which can be integrated to give the result quoted in section [2T 



Appendix B. Solution to the A^-impurities problem 

In solving the Schrodinger equation for ALEs in the Born-Oppenheimer (BO) 
approximation, we must treat the problem of N spin-polarised fermions on a lattice 
(the conduction band atoms), moving in a potential provided by the N holes, which 
essentially act as impurities. The BO Hamiltonian from can be decomposed into 
a sum of single-atom Hamiltonians /i„[R]. As atom-hole interaction is on-site, we can 
write 

N ^ 

^soE] = $^/in[R], /i„[R] = -JpAr„ -C/^ (5r„R„, 
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Each of the single-atom Hamiltonians /i[R] has eigenstates Pfe(r;R), obtained from 
the solving the eigenvalue equations 

/i[R]pfc(r;R) = £;fc[R]pfc(r;R) (B.2) 

where k runs over all bound and scattering solutions. In general, depending on 
dimensionality and the ratio of U/Jp, there will be up to N bound state solutions 
for N holes for any configuration R. If the impurities approach each other too closely 
(on a length scale dependent on [//| Jp|), some of the bound states will disappear. In 
the following, we assume that the system is sufficiently dilute, that N bound states 
will always exist. For N = 2 this will always be the case provided U/\Jp\ > 2. 

From the /9fc(r;R)'s, the exact many-atom wavefunctions 0c((r;R), can be 
constructed by forming the Slater determinant: 

Pfei(ri;R) •■• Pfei(rjv;R) 



.(r;R) 



(B.3) 



Pfe«(ri;R) •■• Pfeiv(rAr;R) 

The wavefunction (r ; R) is antisymmetric in r and solves (fT5|) . The index a thus 
is a multiindex, with a = (fci, . . . , fcjv) and the condition ka ki,, for a ^ h, and the 
BO-energies Ea-BO of (|15p are obtained from summing the single-atom energies: 

{/£l,...,fclv} = Q 

We can then obtain the explicit form of the single-atom wavefunctions pfe(r;R), 
i.e. the solutions to (jB.2p . As the single atom Green's function G(r',r;R) has the 
spectral decomposition G(r',r;R) — p^(r'; R)/9fc(r; R)/(i? — i?fe[R] + irj), (we 
do not write the i?-dependence of the Green's function explicitly in the following), 
knowing G(R„, r; R) gives access to both eigenenergies Ek [R] and wavefunctions (from 
poles and residues respectively). 

As the full single atom wavefunction obeys 

G(r', r; R) = Go(r', r) -f V Go(r', r')V{r')G{r' , r; R) (B.5) 



r' 



where ^(r') = —U J2n=i ^r'TL,, 

r(v'v\-rP[ d^fc exp [ik(r' - r)] 

The set of G(R„, r; R)'s can be determined from the linear system of equations 

A[i?,R]G(r;R) = Go(r), (B.7) 

where 

{A[E,m)nrn = ( ^ ^^."i^^' ^ T " (B.8) 

[ C/Go(Rri,R™) TO /= n 

and G(r;R) = (G(Ri, r; R), . . . , G(RAr, r); R), Go(r) = (Go(Ri, r), . . . , Go(Rjv, r)). 
The eigenenergies are determined from the condition 

det(A[i;fc[R],R]) = (B.9) 

and we assume here that this has N bound-state solutions. We thus see immediately 
that any TV-hole single-atom bound state wavefunction can always be written as a 
linear combination of functions that have the form of single-hole bound states. 
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For example, for = 2 (|B.7p is solved by 

(1 + ^Go(Ri,Ri))Go(Ri,r) - [/Go(Ri, R2)Go(R2, r) 
^^""^'"'-^ (l + C/Go(0,0))^-C/2G^(Ri,R2) 



_ (l + t/Go(R2,R2))Go(R2,r)-[/Go(Ri,R2)Go(Ri,r) 
^^""^'"'-^ (l + C/Go(0,0))^-C/2Gg(Ri,R2) 

where the eigenenergies i?±[Ri,R2] are implicitly determined from the denominator 
being zero, 

(1 + ;7Go(0, Q)f = C/2g2(Ri, R2), (B.12) 

and the residue yields two wavcfunctions, 

p±(r;R)-G(Go(i?i,r)±Go(i?2,r)), (B.13) 

where G carries the normalization. 

The lowest energy BO-potential is given by i?+[Ri, R2] + [Ri, R2], examples 
of which are plotted in figure |4] for the ID case. 



Appendix C. Description of a many-exciton condensate 



The condensed state of the excitons, in which a macroscopic number of them 
occupies the single-exciton groundstate, breaks the symmetry of the Hamiltonian 
associated with number conservation. Macroscopic occupation of the groundstate 
= Ao,o(k, — k)cj^dt,^|0) imphes that (cj^rfl,^) 7^ 0. A way to introduce 
symmetry breaking into the Hamiltonian is by applying a Bogoliubov canonical 
transformation to the particle and hole operators [TJ [5] : 



Gk = r>Ck27^ = UkCk + I'krfLk (C.l) 

Dk = Prfkl?^ ^ ukdk - WkcLk (C.2) 

with Wk+^k = 1 to maintain anti-commutation relations. Here Gk, i?k correspond 
to new quasiparticle and quasihole operators, that have vanishing occupation number 
in the condensate groundstate of the system, 

(G^Gk) - {dId^) = 0, (C.3) 

which automatically imposes {c\a'^_-^^ ^ 0. Inserting the inverse of (jC.ip into and 
rearranging to obtain a normal-ordered form again, the Hamiltonian decomposes into 
a sum of three terms: 

H — Hconst + H2 + Hint (C4) 

These terms have the following structure: 



- Hconst is a constant containing no operators 



- H2 is given by 

H2=Y1 4'^^k^^k + 4'^^k^k - FuiClol^ + D^^C^) (C.5) 

k 
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with 



1 

-Zwkl'k E ^l-k"l^^l 

1 

= 44 - K-l + E [4 (^0 + ^o" - ^1^) + vl {VU -Vo- V^f )] vt 
1 

-2ukl'k E ^i-k^i^i 



(C.6) 



(C.7) 



K + 4 - E (2^0 + ^o" + ^o' - - VU) v. 
1 

+ ("k - ^'k) E ^l-ktil""! 



(C.8) 



The last term in (jC.Sp corresponds to a process where electron-hole pairs with 
total momentum equal to zero are being spontaneously created from and annihilated 
into the condensate. 

- Hint is the transformed interaction part of the Hamiltonian and contains all 
possible quartic combinations of particle and hole operators. Among other processes, 
these correspond to the creation and annihilation of two pairs of quasiparticles and 
quasiholes with total momentum zero. 

H re-expressed in the new operators can thus be unstable with respect to the 
spontaneous creation and annihilation of free quasiparticle-quasihole pairs. To rectify 
this, Keldysh and Kozlov demand that Uk and vy^ be chosen such that they satisfy the 
constraint Ij 

(4^-k> = (Cki?-k) = 0. (C.9) 
The average here is performed over the exact groundstate, and thus Uk and Uk need 
to be determined from the full transformed Hamiltonian, subject to (jC.9p . which is a 
difficult task. 

It is simpler to satisfy constraint (|C.9p just for the groundstate of H2 and disregard 
Hint for the moment. This is equivalent to demanding that Fk = for all k [2], 
which in turn is equivalent to Uk and Vk satisfying mean-field (j26p - (|27p . In this 
approximation, only particle-hole pairs at opposite momenta are correlated. The 
Bogoliubov transformation that generates these correlations is 



V 



exp 



-^A(k)(4dt_k-rf-kCk) 

k 



(C.IO) 



where Wk = cos{y/nexA(k)), Vk = sm(y/n^ A{\i)) . When applied to the vacuum 
particle-hole vacuum |0), V generates the coherent state to the operators cj^dlk' '^-kCk 
which automatically satisfies (|C.3[) and (|C.9p . The new groundstate is consequently 
given by 

It/-) = v\o) - []K + t'k4dLk)lo), (c.ii) 

k 

of the same form as the BCS-groundstate. The dispersion relations of the single- 
quasiparticle excitations above this groundstate are given by (|C.6p and (|C.7p . 
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